Competing mechanisms for singlet— triplet transition in artificial molecules 
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We study the magnetic field induced singlet /triplet transition for two electrons in vertically- 
coupled quantum dots by exact diagonalization of the Coulomb interaction. We identify the different 
mechanisms occurring in the transition, involving either in-plane correlations or localization in 
opposite dots, depending on the field direction. Therefore, both spin and orbital degrees of freedom 
can be manipulated by field strength and direction. The phase diagram of realistic devices is 
determined. 
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Atomic-like phenomenology, ensuing from the discrete 
density of states, has been predicted and demonstrated in 
semiconductor quantum dots (QDs), such as shell struc- 
ture P, @, fine structure due to exchange interaction 
(Hund's rule) 0, and Kondo physics Q; hence QDs are 
often termed artificial atoms. Carriers can be injected 
one by one into the system in single-electron transport Q 
or capacitance Q experiments, based on the Coulomb 
blockade || phenomenon, and the energy required to add 
one electron can be measured if the electrostatic screen- 
ing is poor and the thermal smearing is low. 

Coupled QDs extend to the molecular realm the simi- 
larity between natural and artificial atoms [g, Q; here, 
inter-dot tunnelling introduces an energy scale which 
may be comparable to other energy scales of the sys- 
tem, namely, single-particle confinement energy, carrier- 
carrier interaction, and magnetic energy. In contrast to 
natural molecules, where inter- nuclear coupling is fixed 
by the balance between nuclear repulsion and electro- 
static attraction mediated by valence electrons, in such 
artificial molecules (AM) all energy scales, including 
inter-dot coupling, as well as the charging state of the 
system can be controlled to a very high degree by device 
engineering and/or external fields §. 

A typical AM consists of a disc-like region obtained 
from coupled two-dimensional quantum systems, such as 
two quantum wells (vertically coupled QDs). As in sin- 
gle QDs, electronic states can be easily manipulated by a 
magnetic field B±, perpendicular to the plane of the QDs, 
which drives the system from a low-correlation (low- 
field) regime to a highly correlated (high-field) regime 
by changing the single-particle splittings The study 
of electronic states of few electrons in AMs 0, 0] has 
become a topic of increasing interest, partially due to 
possible implications for the implementation of scalable 
solid-state quantum gates, with the quantum bit of infor- 
mation coded either in the electron charge or s P m E3 
degree of freedom (DOF). 

It should be noted that in AMs carriers are not only 
electrostatically coupled, but also have their spin inter- 
laced when tunnelling is allowed , since electrons with 



opposite spin may tunnel into the same dot if the intra- 
dot Coulomb interaction is not too large; the same pro- 
cess is obviously prohibited for electrons with parallel 
spins. This two-electron dynamics may be described by 
an effective Heisenberg Hamiltonian H = J(B)si • S2 be- 
tween spins Si and S2 [14). with singlet and triplet config- 
urations separated by a field-dependent exchange- energy 
gap J(B) = E t — E s , which is positive at zero field 
One convenient way to control inter-dot tunnelling, and, 
hence, effective spin-spin interaction J, is by applying a 
magnetic field with a finite component in the plane of 
the QDs, i.e., perpendicular to the tunnelling direction 
Controlling tunnelling by Bu has the ad- 
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vantage that other energy scales and, in particular, the 
Coulomb interaction are practically unaffected. However, 
few studies are devoted to this field configuration, which 
lacks the cylindrical symmetry which can be exploited 
in the vertical field arrangement. On the other hand, 
controlling J(B) in AM is crucial for the proposed im- 
plementation of scalable quantum gates |l4j] . 

In this paper we study the exchange energy for two 
electrons confined in AMs in a magnetic field of arbi- 
trary direction. This is performed by a fully numerical, 
real-space approach which allows to account for the com- 
plexity of realistic samples; the carrier-carrier Coulomb 
Hamiltonian is diagonalized exactly within a large single- 
particle basis. We show that the field drives the system 
from an uncorrelated regime, where the singlet state is 
stable, to a strongly correlated one, where triplet order- 
ing is favored; however, the transition occurs by different 
mechanisms, whether the field is in the vertical or in the 
in-plane direction. 

We consider two electrons in a general QD structure. 
Carriers are described by the effective-mass Hamiltonian 
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FIG. 1: Single-particle energy levels for a GaAs AM in a 
magnetic field. Solid and dotted lines represent the FD states 
Enm induced by a strictly vertical field for S and AS levels, 
respectively, s and p shells (see Ref. l2lT) are indicated. Dots 
represent calculated energy levels for a total field of 8 T, ro- 
tated from 0° to 40° with respect to the AM vertical axis. 
Sample parameters are as follows: Lw = 10 nm, Lb = 3nm, 
Vo = 300 meV, and hula = 10 meV. 



with N = 2. Here to*, e* , and g* are the effective 
mass, dielectric constant, and g— factor, respectively |l7| . 
Equation neglects non-parabolicity effects, but oth- 
erwise includes the full 3D nature of the quantum states 
in realistic samples, such as layer width and finite band 
offsets, by the effective potential V(r). Our numerical 
approach consists in mapping the single-particle terms 
in a real-space grid, leading to a large sparse matrix 
which is diagonalized by the Lanczos method. Single- 
particle spin-orbitals are then used to build a basis of 
Slater determinants for the A^-particle problem, which 
is then used to represent the two-body term, in the fa- 
miliar Configuration Interaction approach Q. Coulomb 
matrix elements are calculated numerically. The ensuing 
matrix, which can be very large, is again sparse and can 
be diagonalized via the Lanczos method as well 0] . 

In the following the potential V(r) describes two iden- 
tical vertically coupled disk-like QDs. As usual for this 
type of samples which have very different confinement 
energies in the growth and in-plane directions, we sep- 
arate the potential as V(x,y) + V(z), where V(z) rep- 
resent two symmetric quantum wells of width Lw sepa- 
rated by a barrier Lb and conduction band mismatch Vq. 
We perform the common choice of a parabolic in-plane 
confinement {l/2)m*u!n(x 2 + V 2 ). as this has proved to 
be quantitatively accurate Note, however, that our 
numerical approach does not assume any symmetry; in 
particular, the vector potential A(r) is not limited to 
describe z-directed field. 

In a QD with parabolic in-plane confinement and 
strictly perpendicular magnetic field, single-particle 



states are given by the Fock-Darwin (FD) states (see, 
e.g., Ref. 1!)}. with energies e nrn = htt(2n + \m\ + 
1) — (hu> c /2)m, n and m being the principal and az- 
imuthal quantum numbers, respectively. The oscillator 
frequency is fl — y^oj 2 + ojf/i, with the cyclotron fre- 
quency tu c = eB/m*c. In symmetric AMs we have two 
such ladders of energy levels, associated with the sym- 
metric (S) and anti-symmetric (AS) states arising from 
the double-well potential in the growth direction, rigidly 
separated by a splitting Asas (see Fig. EJ. 

We next consider the effect of a magnetic field with a 
finite in-plane component Bn . As shown in Fig. ^ when 
the angle 9 between a fixed |B| and the z axis is increased, 
the energy levels no longer correspond to the FD states 
at the corresponding B±. Indeed, the splitting between 
S and AS levels decreases with increasing 9 |2fJ, which 
shows that an in-plane component of the field suppresses 
the tunnelling; note that this effect is larger for higher 
levels. It is important to stress that the in-plane field 
can meaningfully affect the motion along the growth di- 
rection if LUc — &Bn/m*c ~ AgAs- Similar effects are 
much harder to achieve in single QDs, due to the large 
single-particle gaps induced by the single quantum well 
confinement. 

As discussed in more detail in the following, the reduc- 
tion of the energy gap between the s and the p shells j2]| 
(see Fig. ^| strongly reduces the single-particle energy 
of the triplet state with respect to that of the singlet: 
the perpendicular field thus promotes the singlet-triplet 
crossing. This transition results in an enhancement of 
the in-plane correlation of the two-electron groundstate 
and in the spin-polarization of the system Il9ll22| . arising 
from the exchange (orbital) interaction |23j. Note that 
this mechanism only involves the in-plane DOFs, and is 
therefore present in both single and coupled QDs; in or- 
der to observe some marked differences in the behavior 
of the two systems, one needs to excite the motion along 
the growth direction z. 

Figure Efa) shows the single-particle levels as a func- 
tion of the in-plane field Bh. The energy levels come 
in shells with S and AS character, but the degeneracies 
which are present at Bn = are removed by a finite field, 
as the axial symmetry of the system is lost. Therefore, 
the single-particle wavefunctions do not have a well de- 
fined angular momentum, and are now S or AS only with 
respect to a 180° rotation about the axis parallel to B. 
Besides, as the field is increased, the S and AS levels ap- 
proach each other, since the tunnelling is progressively 
suppressed [lfij . 

In Fig.[2Ib) we show the lowest two-particle levels, and 
schematically indicate the main components of the cor- 
responding wavefunctions in terms of S and AS single- 
particle states. At low B\\ the ground- and the first 
excited-state have a singlet and a triplet character, re- 
spectively. As Bn is increased, the energy gap J is sup- 
pressed: indeed, singlet and triplet states have the same 
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FIG. 2: Energy levels vs in-plane field at B± = for the 
same AM of Fig. (a) Single-particle levels, with indication 
of the S/AS character at low field, (b) Two-electron levels. 
Insets: main components of the wavefunctions in terms of S 
(left boxes) and AS (right boxes) single-particle states. 

orbital energy, while the Zeeman term favors the latter 
(in the field range of Fig.^the Zeeman contribution can 
hardly be distinguished). As shown in the insets, the 
transition occurs with the maximal, Coulomb-induced 
mixing of the S and AS states, which is favored by the 
vanishing of Asas at large fields. In other words, in- 
creasing Bu the singlet state evolves from a nearly pure 
S state to a fully entangled state in the S/AS basis. Note 
also that, contrary to the one occurring at large B±, here 
the transition is associated to the correlation along the 
growth direction, i.e., with the two electrons sitting on 
opposite QDs, as we will show below. 

The ability to control both the exchange energy J and 
the effective Hilbert-space structure is indeed pivotal to 
the QD-based implementations of quantum-information 
processing Q; besides, J is of direct experimental in- 
terest, for it can be probed by single-electron excitation 
spectroscopy [24|. In Fig. EI a) we show the calculated 
exchange energy as a function of the in-plane field at 
different values of B±_ and for a weaker parabolic con- 
finement (Ag^s < ?koo). The positive/negative J region 
is the stability region for singlet/triplet states. Figure 

a) shows that an increase in B±_ monotonically (i) re- 
duces the singlet stability range with respect to Bn , and 



FIG. 3: (a) Exchange energy J vs in-plane field at selected 
vertical fields, for a GaAs AM. Sample parameters are as in 
Fig-Q but for a weaker lateral confinement hcoo — 4meV. For 
clarity we show the best fitting curves from a large number of 
calculated points. Numerical inaccuracies may result in ±0.3 
meV shift from the curves only for the highest Bj_. (b) Cal- 
culated singlet /triplet phase diagram. The line is a guide to 
the eye through the calculated points [26| . Insets show the 
singlet (solid line) and triplet (dashed line) conditional prob- 
ability near the transitions, defined as \tp(po, zo, p\ z)\ 2 ; tjj is 
the two-electron wavefunction, with p the in-plane coordinate 
with respect to the vertical axis of the cylindrical QD. The 
reference electron (black dot) is fixed at zo = —7.5 nm, at an 
in-plane position \po\ — 4.4 nm; the conditional probability 
is then plotted along an axis parallel to z and crossing the 
QD plane at a position diametrically opposed to the refer- 
ence electron. Left inset: B± — 4T, By = 0T. Right inset: 
B ± = 0T, By =9T. 

(ii) enhances the ferromagnetic (J < 0) behavior in the 
considered range of Bu values. These features are sum- 
marized in the phase diagram shown in Fig. |3{b). 

A closer inspection into our results shows that different 
mechanisms are involved in the singlet /triplet transition, 
depending on the field direction. At zero field the singlet 
state mainly corresponds to both electrons occupying the 
(s,S) orbital, while a minor contribution from the (s,AS) 
orbital gives rise to the spatial correlation in the z direc- 
tion. On the contrary, all the dominant configurations in 
the triplet state involve S states (see the conditional prob- 
ability in Fig. Efb)). B± leaves unaffected the z DOF, 
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while it energetically lowers the p (to = 1) state with re- 
spect to the s one. The positive single-particle contribu- 
tion to J is therefore reduced, until it is compensated by 
the negative contribution arising from the Coulomb en- 
ergy. To summarize, the singlet/triplet crossing induced 
by B± is mainly connected with the in-plane dynamics, 
while it leaves unaffected the motion in the growth direc- 
tion and the double-occupancy probability of each dot. 

The main effect of B\\ (right inset), instead, is that of 
suppressing the energy splitting resulting from the inter- 
dot tunnelling. This clearly favors the occupation of the 
AS states, and therefore vertical correlations for both 
the singlet and the triplet states set in; in both cases, 
the two electrons tend to localize in opposite dots, and 
the importance of the spins' relative orientation vanishes 
with the double occupancy probability. Indeed the ex- 
ponential vanishing of J represents the clear fingerprint 
of the regime where the double-occupancy probability is 
suppressed. It should be noted that this is not a single- 
particle effect, since it does not imply, nor require, the 
complete suppression of the tunnelling. 

The results reported in Fig. [21 show that these two dif- 
ferent mechanisms interfere with each other in a non- 
trivial manner. The presence of the perpendicular com- 
ponent B± favors the single-triplet crossing and the fer- 
romagnetic phase, while it opposes the suppression of 
the double occupancy and the resulting singlet-triplet 
degeneracy (apart from the Zeeman term). Such inter- 
play arises from the 3D nature of quantum states in the 
AM: in fact, in the considered range of physical parame- 
ters (Tiujo ~ Asas), the magnetic field can strongly affect 
both the in-plane (intra-dot) and the vertical (inter-dot) 
DOF. 

An adiabatic manipulation of J by means of magnetic 
(and electric) fields has been proposed in order to imple- 
ment the two-qubit gates in electron-spin based quantum 
computers |l4j . The rest condition within such scheme 
would correspond to the suppression of J and of the over- 
lap between electrons localized in adjacent QDs, where 
both conditions should be induced by a static magnetic 
field. In this perspective, our findings suggest that (i) 
the -By (rather than B±) component of the field and the 
exponential suppression (rather than the crossing point 
from J > to J < 0) are required; (ii) the presence 
of a field component perpendicular to the static one (as 
required, e.g., for the single-spin rotations) should be si- 
multaneously taken into account in order to determine 
the suited range of physical parameters. 

To summarize, we have theoretically investigated the 
dependence of the singlet and triplet states of two elec- 
trons in AMs on external magnetic fields of arbitrary 
direction. Our computational approach allows to fully 
account for the different physical mechanisms underlying 
the singlet /triplet transitions which are due to the paral- 
lel and perpendicular components of the field, as well as 
for the non-trivial interplay between the vertical and the 



in-plane correlation effects that they induce. The perpen- 
dicular component of the field does indeed facilitate the 
transition from the anti-ferromagnetic to the ferromag- 
netic phase which is induced by the parallel component, 
but at the same time it opposes the carrier localization 
and correlation properties that the latter tends to induce 
in the AM. 
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